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1.0  INTRODUCTION 

Plume  diagnostic  measurements  conducted  at  Arnold  Engineering  Development  Center 
(A£DC)  on  several  solid-propellant  rocket  exhausts  have  revealed  phenomena  which  may 
provide  insight  into  the  characterization  of  these  plumes.  Two  of  these  phenomena,  aluminum 
chloride  (AlCl)  absorption  and  searchlight  emission,  are  highlighted  in  this  report.  Modeling 
of  each  phenomenon  and  comparison  of  predictions  with  plume  measurements  acquired  during 
simulated  altitude  chamber  tests  are  also  presented. 

A  discussion  is  presented  in  Section  2.0  on  the  analysis  of  the  AlCl  absorption  feature 
observed  in  aluminized  solid-propellant  rocket  plumes.  The  results  of  a  line-by-line  calculation 
are  discussed  and  compared  with  high-resolution  ultraviolet  (UV)  spectral  data. 

In  Section  3.0,  selected  plume  radiation  measurements  are  reviewed,  and  evidence  of  search¬ 
light  emission  is  presented.  Two  modeling  approaches  are  discussed  and  compared  with  motor 
test  data. 


2.0  ALUMINUM  CHLORIDE  ABSORPTION  ANALYSIS 

2.1  INTRODUCTION 

An  interesting  feature  in  the  mid-ultraviolet  spectrum  of  the  radiation  from  aluminized 
solid-propellant  rocket  exhausts  has  been  noted  several  times  in  unpublished  reports.  This 
feature  is  a  deep  absorption  near  262  nm  of  the  otherwise  continuum  spectrum  typical  of 
hot  aluminum  oxide  particles.  It  has  been  attributed  to  the  AlCl  A-X  transition,  although 
the  spectral  resolution  was  not  sufficient  to  resolve  bandheads  (Ref.  1).  Recently,  spectra 
were  obtained  having  sufficient  resolution  to  identify  the  vibrational  bands,  and  these  will 
be  reported  here.  Curiously,  spectra  obtained  from  an  onboard  spectrometer  viewing  an  in¬ 
flight  rocket  plume  having  an  aluminized  propellant  also  contained  this  feature  (Ref.  2). 

One  obvious  question  arising  from  such  an  inherent  spectral  feature  is  whether  it  can 
be  used  to  diagnose  any  properties  of  the  rocket  exhaust  flow.  To  that  end,  a  modeling  effort 
was  undertaken,  and  this  paper  reports  on  the  success  of  that  work. 

2.2  SPECTRA  DESCRIPTION 

When  a  moderate-resolution  (0.8  nm)  ultraviolet  scanning  monochromator  is  pointed 
normal  to  the  axis  of  an  exhaust  plume  from  an  aluminized  solid-propellant  rocket  motor 
being  fired  in  an  altitude  test  cell  (Fig.  1),  the  spectra  shown  in  Fig.  2  result.  The  continuum 
on  either  side  of  the  feature  at  262  nm  and  the  thermal  behavior  of  the  radiation  in  the  visible 
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and  infrared  regions  suggests  that  the  feature  is  caused  by  the  absorption  of  a  gaseous  specie. 
A  quick  examination  of  the  predicted  products  of  combustion  of  ammonium  perchlorate 
with  aluminum  and  hydrocarbon  solids,  along  with  a  trip  through  Pearse  and  Gaydon’s  (Ref. 
3)  atlas  of  spectra,  reveals  that  the  most  likely  absorber  is  AlCl.  To  investigate  this  identification 
more  thoroughly,  higher-resolution  spectra  were  obtained  during  a  subsequent  test.  Because 
of  the  small  vibrational  (approximately  450  cm~  0  and  rotational  (approximately  0.25  cm~  *) 
constants  of  both  AlCl  states,  a  spectral  resolution  of  about  0.1  nm  was  required  to  produce 
the  band  structure.  The  resulting  spectrum  (Fig.  3)  indeed  does  show  the  expected  structure, 
and  many  of  the  bands  listed  in  Pearse  and  Gaydon  (Table  1)  are  identified. 

2.3  AlCl  MOLECULE 

Aluminum  chloride  sublimates  at  temperatures  above  460  K.  A  calculation  of  the  plume 
properties  using  the  JANNAF  standard  plume  flow-field  codes  (Ref.  4)  gives  the  solid  and 
gaseous  tempierature  profiles  of  Fig.  4a  for  this  particular  rocket  plume.  The  minimum 
predicted  temperature  is  about  600  K  so  that  the  AlCl  is  gaseous.  Furthermore,  the  predicted 
mole  fraction  of  AlCl  is  about  2.6  percent  and  should  follow  the  static  pressure  radial  profile 
(Fig.  4b). 

The  molecular  constants  of  the  two  electronic  states  (A^n  and  X^),  measured  by  several 
investigators  and  reported  critically  by  Huber  and  Herzberg  (Ref.  5),  are  given  in  Table  2. 
The  values  are  for  Al^^Cl,  although  Al^^CI  also  naturally  occurs.  These  two  isotopes  occur 
in  a  ratio  of  75.5  to  24.5.  The  relatively  small  vibrational  constant  of  the  ground  state  (480 
cm~‘)  leads  to  a  vibrational-rotational  fundamental  at  a  wavelength  of  about  21  lan. 
(Observation  of  this  LWIR  emission  band  from  rocket  plumes  has  not  been  reported.)  The 
constants  of  the  state  have  been  re-measured  recently  in  the  infrared  bands,  and 
improvements  in  accuracy  have  been  reported.  These  data  will  be  introduced  shortly.  The 
rotational  constants  for  this  closely  packed  molecule  are  also  quite  small  (approximately  0.254 
cm~^),  and  the  overlapping  branches  make  resolution  of  the  lines  nearly  impossible.  The 
vibrational  correction  to  the  rotational  constant  given  as  footnotes  in  Ref.  5  is  a  necessity 
for  proper  modeling  of  this  band.  Also,  because  of  the  small  values  of  Be,  rotational 
quantum  numbers  out  to  at  least  100  must  be  used  at  the  temperature  of  interest. 

2.4  RADIATIVE  TRANSFER  MODEL 

The  preliminary  model  adopted  for  this  study  is  illustrated  in  Fig.  5.  It  is  supposed  that 
the  core  of  solid  particulates  provides  a  continuum  source  for  transmission  through  a  layer 
of  gas  containing  the  absorbing  AlCl.  The  continuum  radiation  from  the  particulates  was 
curve  fit  by  joining  the  spectra  at  both  extremes  of  the  absorption  feature,  and  served  as 
a  source  for  the  absorption  calculations.  The  gas  layer  was  considered  uniform  for  these 
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preliminary  calculations,  with  temperature  and  number  density  as  parameters.  The  absorption 
spectra  were  modeled  as  Voigt  lines  with  the  spectrometer  triangular  slit  function  considered 
to  be  much  greater  than  a  line  width, 

A  complete  solution  to  the  monochromatic  radiative  transfer  equation  is  desired.  The 
equation  is  written 


^  =  Se  Ntot  L-  -  Sa  Ntot  K  I,  (1) 

dx 


where 

Ntot  =  Total  number  density  of  emitter/absorber 

L,  =  Line  shape  function 

Se  =  Molecular  emission  strength 

Sa  =  Molecular  absorption  strength 

I,  =  Spectral  radiance  (i.e.,  W/sr-cm^  -  cm“*)  at  wavenumber  v 

X  =  Spatial  coordinate  within  the  media 

The  boundary  condition  at  x  -  0  is  given  in  this  case  by  the  tit  to  the  continuum  spectra, 
neglecting  the  AlCl  absorption  feature;  i.e.,  the  continuum  background. 

Conceptually,  the  calculation  involves  3  steps: 

1.  Determine  the  positions  of  all  spectral  lines.  Determine  the  AlCl  emission  and 
absorption  strengths  for  each  line  based  on  a  rotational  and  vibrational 
temperature.  This  yields  a  spectral  line  atlas  of  line  strengths  versus  line  center 
wavenumber  for  aU  lines. 

2.  Given  Sa  and  Se,  and  Njoti  solve  Eq.  (1).  Here,  only  one  spatial  zone  is 
considered  so  that  integration  of  Eq.  (1)  is  simplified. 

3.  Convolve  the  result  of  step  2  for  each  line  with  spectrometer  response  function. 

In  this  case  a  triangular  slit  function  suffices  since  the  actual  line  width  is  much 
smaller  than  the  instrument  monochromatic  response  function. 
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2.4.1  Part  1 


The  strength  of  each  emission  line  is  given  by 

hci'  .  Sj'j' ' 


Se  =■ 


4n 


2J'  +  1 


x(A;v';J') 


(2) 


where 

Av'v' 

Sj’j" 


=  Einstein  coefficients,  (Ref.  7) 


=  Honl-London  factor  (Ref.  8) 


=  Wavenumber 


x:(A;v';J')  =  Fractional  population  of  the  A*n  state  for  v'  and  J'  level,  normalized 

by  the  partition  function,  i.e., 


x(A,v',J')  = 


gAv'j  expl-  -^E(A,v',J') 


Eacxp|--^E,. 


(3) 


where  gi  is  the  degeneracy  factor  of  the  given  state.  The  absorption  strength  is  similarly  given 
by 

Sa  =  hj/Bv'v"  •~f^{x(X.v".J")  -  x(A.v'.J')}  (4) 


with  the  Einstein  B  coefficient  as 


Bv 


1 


Srho'^ 


(5) 


The  band  strength  (or  lifetime)  of  the  AlCl  A^II  v  ~  0  state  has  been  measured  (Ref, 
6)  to  be  6.4  x  10"^  sec,  where  no  effort  was  made  to  separate  the  two  isotopic  states.  The 
lifetime  was  calculated  to  be  S.17  x  10~^  sec  by  Langhoff  and  Bauschlicher  (Ref.  7),  who 
also  calculated  relative  band  strengths  from  first  principles.  The  transition  probabilities  for 
each  significant  vibrational  transition  are  given  as  Einstein  coefficients  in  Table  3,  where 
the  calculated  lifetime  was  used.  Note  the  relatively  strong  values  for  the  Av  =  0  sequence. 
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indicating  that  these  will  be  the  most  evident  in  the  structure,  as  also  reported  in  the 
experimental  analyses  (Ref.  3).  These  calculations  of  band  transition  probabilities  were  made 
for  Al^^Cl,  with  the  assumption  that  the  results  are  the  same  for  Al^’Cl. 

Energies  of  individual  states  are  given  by 

Ei  =  Tj  =  Evj  +  Er;  (6) 

where  Ti  is  electronic  state  energy,  Evj  is  the  vibrational  state  energy,  and  Erj  is  the  rotational 
energy  of  quantum.  Now,  in  terms  of  the  vibrationai  quantum  number  v,  the  rotational 
quantum  number  J,  and  the  electronic  angular  momentum  quantum  number  A, 

Evj  =  cue  (v  +  Vi)  -  Wg  x*  (v  +  •/i)2  +  Ug  y,  (v  +  Yi)^  (7) 

Er,  =  Bv  [J(J  +  1)  -  A2]  -  Dv[(J(J  +  1)  -  A2)2]  (8) 


with 


By  =  Be  Qe  (V  +  ^i)  +  a't  (v  +  !4)2  +  a'^'  (v  +  1/2)3 


Dy  =  Dg  +  i8e  (v  +  Yi) 

For  mass  AP^ci,  all  constants  are  modified  by  the  ratio  g  =  y/iSTfT 

Bf  ^  =  Bf  >  g2 


af '  =  4''’^  e3  etc. 


The  molecular  constants  for  the  AiJI  state  were  taken  from  Huber  and  Herzberg  (Ref.  S) 
and  are  listed  in  Table  2.  For  the  ground  electronic  state,  Bernath’s  measurements  use  the 
Dunham  expression  (Ref.  8)  for  the  sum  of  the  vibrational  and  rotational  energy, 

Tyj  =  E  Y,„  (v  +  Yiy  [J(J  +  l)r 

i,m 

The  values  of  Yim  were  given  to  the  authors  by  Bemath  and  appear  in  Table  4.  With  the 
transition  energies  defined,  line  positions  are  obtained  by  subtraction 


V  =  E(A,vM')  -  E(X.v",J") 


A  line  atlas  is  built  up  from  these  expression  for  v',  v"  out  to  8  and  J' ,  J"  out  to  200 
for  both  isotopes. 
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2.4.2  Part  2 

The  line-by-Iine  model  then  solves  Eq.  (1)  for  each  line  in  the  atlas.  A  Voigt  profile  for 
each  line  is  assumed.  The  Sb  and  Sa  are  calculated  at  a  minimum  of  ten  spectral  frequencies 
between  spectral  line  centers  by  summing  contributions  of  all  relevant  lines  at  that  position. 
Here  v'  ,v"  out  to  8  and  J'  J"  out  to  200  are  used  to  generate  a  large  number  of  stored 
numbers.  Thus,  for  20,000  spectral  lines,  at  least  200,000  radiances  and  transmittances  are 
generated.  The  program  also  does  spatial  integration,  although  for  this  case  (1  zone)  the 
integration  is  trivial. 

2.4.3  Fail  3 

The  spectrometer  is  modeled  by  a  triangle  function  of  2  A  full  width  at  half  maximum. 
The  convolution  of  the  integral 

H(.»)  =  i  (10) 

is  done  numerically  using  the  trapezoid  rule,  where  H(i')  represents  the  instrument  reading 
at  V  and  G(r')  represents  the  instrument  function, 

2.5  MODELING  RESULTS 

The  approach  taken  was  to  vary  the  vibrational  and  rotational  temperatures,  the  column 
density,  and  the  column  length  in  Eq.  (2)  until  the  closest  fit  to  the  experimental  data  was 
achieved.  The  temperature  and  pressure  were  varied  about  the  predicted  averages  along  the 
column,  (see  Fig.  4)  and  the  predicted  mole  fraction  of  AlCl  was  used  initially.  It  was  first 
necessary  to  calculate  the  constants  accurately  for  each  vibrational  state  of  each  of  the  isotopic 
molecules.  This  proved  to  be  very  sensitive  in  obtaining  a  fit  to  the  vibrational  bandheads. 
Then,  the  number  of  J  values  necessary  to  fit  the  absorption  level  was  determined. 

Only  two  results  will  be  shown;  First,  the  result  using  the  predicted  temperature  and  AlCl 
density  (Fig.  6)  and  the  final  result,  the  best  fit  (Fig.  7).  As  shown  in  Figs.  4  and  S,  the 
parameters  used  for  the  baseline  calculations  were 

Vibrational  Temperature,  Tv  =  800  K 

Rotational  Temperature,  Tr  =  800  K 

Number  Density,  t)  =  2.36  x  10**  crn’’^ 

Absorber  Path  Length,  f  =  5  cm 
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Quite  clearly  from  Fig.  6,  the  fit  of  the  calculation  to  the  data  is  poor.  The  structure 
is  much  smoother  than  the  data,  although  the  same  spectral  resolution  was  used  in  the  model 
as  for  the  spectrometer.  In  the  process  of  iterating  on  temperature  and  density,  a  better  fit 
could  be  found  in  the  central  portion  of  the  absorption,  but  the  structure  in  the  wings  became 
poorer.  It  was  this  lack  of  spectral  fit  that  led  to  use  of  the  more  accurate  data  for  the 
state.  Finally,  since  the  molecular  data  for  the  upper  state  are  somewhat  dated  and  some 
uncertainty  exists  in  the  relative  values  of  the  rotational  constant  for  the  two  isotopes,  the 
vibrational  and  rotational  constants  of  the  A^II  state  were  also  iterated.  Since  this  process 
could  go  on  and  on,  a  final  fit  (Fig.  7)  was  produced.  The  final  parameters  are  as  follows; 

Vibrational  Temperature,  Ty  =  1,000  K 

Rotational  Temperature,  Tr  =  1,000  K 

Number  Density,  rj  =  2.36  x  10‘®  cm"^ 

Absorber  Path  Length,  f  =  5  cm 

and  the  molecular  constants  were  also  shifted: 

«e  450  —  430  cm“* 

cogXc  4.3  —  4.1  cra~> 

Be  0.259  -  0.249  cm"' 

Here,  the  spectral  match  is  fair  and  the  overall  absorption  matches  pretty  well,  so  that  continued 
iteration  appears  unwarranted. 

The  final  result  is  reasonable.  Because  of  the  high  particle  temperature,  it  is  not  surprising 
to  have  a  gas  temperature  somewhat  larger  than  predicted  (1,000  K  versus  800  K  predicted). 
Also,  because  of  small  particles  in  the  outer  gas  flow,  it  is  not  surprising  that  the  radiation 
from  the  gas  layer  masks  the  absorption  at  the  deepest  absorption  frequency.  One  might 
expect  that  an  arbitrary  shift  of  the  wavelength  for  the  experimental  data  would  produce 
a  better  fit.  However,  attempts  to  make  this  match  were  only  marginally  successful.  The 
simple  physical  model  of  Fig.  5  cannot  be  correct,  since  the  radiation  is  obviously  extended 
and  the  AlCl  exists  all  through  the  plume.  However,  the  modeling  brought  assurance  that 
the  absorber  was  AlCl,  and  only  a  slight  variation  of  temperature  brought  the  model  results 
very  near  the  measured  data. 
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2.6  CONCLUSIONS 

Clearly,  if  the  AlCl  feature  is  to  be  used  for  diagnostics,  laboratory  work  using  a  heated 
absorption  cell  will  be  required  to  reach  a  better  fit  between  calculated  and  measured  spectra. 
In  a  heated  cell  experiment  the  temperature  and  pressure  parameters  would  be  well  known, 
leaving  only  the  spectrometer  wavelength  calibration  and  the  molecular  constants  unknown. 
The  issues  raised  above  can  then  be  resolved. 

3.0  SEARCHLIGHT  EMISSION  ANALYSIS 
3.1  INTRODUCTION 

Prediction  of  plume  radiative  heat  transfer  to  spacecraft  is  important  for  the  design  of 
vehicle  thermal  protection  systems.  Radiative  heating  Is  of  particular  concern  for  aluminized 
solid-propellant  rocket  motors  because  of  the  intense  thermal  continuum  emission  from  their 
exhaust  plumes.  The  continuum  is  largely  due  to  thermal  emission  from  hot  micron-sized 
aluminum  oxide  droplets  and  particles  in  the  exhaust  plume  with  temperatures  in  the  range 
2,000  -  2,500  K.  Because  of  uncertainties  in  AI2O3  size,  temperature,  and  cmissivity,  a  priori 
code  predictions  are  unreliable.  Spacecraft  designers  have  instead  relied  on  empirical  methods 
to  predict  in-flight  heat  transfer.  This  works  well  for  launch  vrfiicles  because  it  is  usually 
possible  in  a  ground  test  to  measure  the  heat  flux  in  the  desired  locations  and  viewing  angles. 

The  situation  is  different  for  space  motors,  however,  because  these  are  test^  in  altitude 
simulation  facilities  with  close  proximity  to  an  exhaust  diffuser  duct.  It  is  almost  never  possible 
to  observe  more  than  a  meter  or  so  of  undisturbed  plume  flow.  Sometimes  only  a  few  inches 
of  pltune  are  visible.  In  addition,  the  diffuser  restricts  the  achievable  viewing  angles  to  near 
90  deg,  broadside  to  the  plume.  The  viewing  angles  important  for  spacecraft  heating,  by 
comparison,  are  nearer  to  30  deg.  One  is,  therefore,  left  with  the  problem  of  reconstructing 
the  radiative  heat  flux  from  this  limited  database.  Bobco  (Ref.  9)  and  Edwards  (Refs.  10-12) 
describe  the  current  industry  standard  method  for  curve  fitting  plume  data  from  altitude 
tests  in  which  the  axial  decay  in  radiosity  is  measured,  to  a  radiative  exchange  view  factor 
model.  The  principal  difficulty  in  applying  this  method  is  the  need  to  extrapolate  the  full 
spatial  and  aspect  angle  dependences  of  the  plume  radiance  based  on  a  limited  number  of 
radiometer  or  IR  camera  measurements,  all  of  which  are  necessarily  limited  to  the  nozzle 
exit  region  and  to  roughly  90-deg  aspect  angle. 

Recent  measurements  of  several  motors  at  the  AEDC  (Ref.  13)  suggest  that  a  significant 
portion  of  the  broadside  exit  plane  radiance  of  space  motors  may  be  due  to  nozzle  searchlight 
emission.  Specifically,  the  broadside  radiance  near  the  nozzle  exit  plane  appears  to  contain 
significant  contributions  from  the  hot  inner  nozzle  wall.  Although  not  dirertly  in  the  line 
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of  sight,  the  intense  nozzle  radiation  is  scattered  into  the  field  of  view  by  AI2O3  droplets 
and  particles  in  the  plume  (Fig.  8).  There  has,  therefore,  been  great  interest  in  understanding 
both  the  strength  of  this  searchlight  emission  and  its  implications  for  in-flight  heat  transfer. 
This  is  not  easy,  as  radiometers  and  IR  cameras  cannot  /Ter  se  distinguish  between  searchlight 
and  direct  plume  emission.  Instead,  the  identification  of  searchlight  emission  must  be  made 
on  the  basis  of  characteristic  temporal,  wavelength,  and  spatial  properties  of  the  radiance 
data.  The  following  sections  indicate  what  some  of  these  characteristic  properties  of  search¬ 
light  emission  are  and  how  they  occur  in  selected  motor  test  data. 

3.2  SEARCHUGHT  MECHANISM 

This  section  describes  the  radiative  transport  properties  in  rocket  exhaust  plumes  and 
shows  why  these  should  naturally  lead  to  a  strong  searchlight  effect.  To  dispel  possible 
confusion,  it  should  first  be  mentioned  that  the  term  “searchlight’*  is  a  misnomer.  It  was 
originally  proposed  on  the  mistaken  assumption  that  the  radiation  in  question  originated 
from  the  extremely  hot  combustion  chamber  (Ref.  14).  Additionally,  the  word  “searchlight” 
implies  a  degree  of  collimation  which  simply  does  not  occur  in  rockets.  Since  it  has  become 
entrenched  in  the  technical  jargon,  however,  we  continue  to  use  it  in  this  report. 

The  searchlight  mechanism  is  a  light-scattering  process  controlled  by  the  radiative  transport 
properties  of  the  plume.  The  scattering  is  due  to  micron-sized  AI2O3  particulates  which 
typically  comprise  30  percent  of  the  mass  flow.  Since  the  particulate  temperatures  bracket 
the  material  melt  point  of  2,327  K,  both  solid  particles  and  liquid  droplets  occur.  Figure 
9  shows  that  these  particles  are  very  efficient  scatterers  of  infrared  radiation.  Specifically, 
this  figure  shows  the  Mie  scattering  albedo  of  a  representative  3-^-radius  liquid  AI2O3 
droplet  between  1-  and  4-^m  wavelength.  It  was  generated  using  the  BHMIE  code  of  Ref. 
15  and  the  real  refractive  index  data  of  Refs,  16-19.  The  spectral  structure  is  due  to  Mie 
resonances  and  would  be  smoothed  out  for  a  distribution  of  particle  sizes.  As  can  be  seen, 
the  albedo  is  close  to  99  percent  over  this  entire  spectral  region.  It  would  be  even  higher 
for  a  solid  particle.  The  directional  properties  of  this  scattering  are  indicated  in  Fig.  10,  which 
shows  the  normalized  scattering  phase  functions  for  the  same  3-^m-radius  droplet  at  1-,  2-, 
and  S-itm  wavelength.  The  scattering  efficiency  is  highest  in  the  forward  direction  (0  deg) 
and  drops  rapidly  with  increasing  angle.  The  inset  indicates  how  this  directionality  affects 
the  intensity  of  searchlight  emission  at  different  aspect  angles.  It  considers  emission  from 
the  nozzle  wall  which  is  scattered  at  a  plume  volume  element  dV.  The  ray  path  labeled  A 
contributes  to  the  broadside  radiance.  Since  the  scattering  angle  is  small,  the  efficiency  is 
high.  Path  B  contributes  to  vehicle  heating.  Since  the  scattering  angle  is  large,  the  efficiency 
is  low.  Based  upon  this  reasoning,  the  searchlight  is  expected  to  contribute  more  strongly 
to  broadside  radiance  than  to  spacecraft  heating. 


15 


AEDC-TR-92-12 


In  addition  to  the  scattering  properties  of  single  particles,  searchlight  is  also  affected  by 
the  aggregate  transport  properties  of  the  total  plume.  The  two  important  parameters  in  this 
regard  are  the  scattering  and  the  absorption  optical  depths,  and  Ta,  rcsp>ectively.  An 
estimate  of  these  parameters  can  be  made  using  the  data  of  Fig.  9,  the  nominal  particle  size 
given  by  the  Hermsen  size  correlation  (Ref.  20),  and  assuming  one>dimensional  source  flow 
as  indicated  in  Refs.  21-23. 

Ji-  =  _!5_  =  (1.38  E  -  4  (11) 

Q.  Qa  2xgRv 


where 

Qa  =  Mie  absorption  efficiency  parameter 
Q,  =  Mie  scattering  efficiency  parameter 

Dt  =  Nozzle  throat  diameter 

m  =  AI2O3  mass  flow  rate 

R  =  Plume  radius 

Q  -  Bulk  density  of  AI2O3 

V  =  Particle  velocity 

All  parameters  are  in  cgs  units.  For  small  to  mid-sized  motors,  the  predicted  scattering  optical 
depth  is  on  the  order  of  unity  or  larger  (i.c.,  optically  thick).  This  means  that  most  photons 
propagated  along  a  diameter  through  the  plume  will  be  scattered  at  least  once  and  possibly 
several  times.  The  absorption  optical  depth,  on  the  other  hand,  is  only  several  percent  (i.e., 
optically  thin).  This  means  that  radiation  can  propagate  through  the  plume  without  appreciable 
absorption.  Another  consequence  of  the  low  absorption  optica!  depth  is  that  the  plume  emission 
will  tend  to  scale  linearly  with  thrust. 

A  simple  gedanken  experiment  further  clarifles  both  the  strength  and  spatial  dependence 
of  the  searchlight.  Consider  a  single  alumina  droplet  suspended  over  a  hot,  flat  plate  with 
unit  emissivity  (Fig.  II).  The  droplet  both  emits  radiation  and  scatters  radiation  from  the 
plate.  One  then  observes  the  total  droplet  radiance  along  a  line  of  sight  paraUel  to  the  surface 


16 


AEDC-TR-92-12 


of  the  plate  and  which  at  no  time  views  the  plate  directly.  If  the  height  of  the  droplet  above 
the  plate  z  is  much  smaller  than  the  size  of  the  plate  L,  the  plate  appears  essentially  infinite. 
The  thermal  emission  from  the  droplet  is 

irriQaB(Td)  (12) 


where 

r  =  Droplet  radius 

Qa  =  Mie  absorption  efficiency  parameter 

B  =  Planck  radiation  function 

Td  =  Droplet  temperature 

The  radiation  emitted  by  the  plate  and  scattered  into  the  line  of  sight  (x  direction)  by  the 
droplet  is  the  searchlight  emission.  It  is  given  by 

irr2  Q,  B(Tp)  I  P(n  -  x)  dO  (13) 


where 

Qs  =  Mie  scattering  efficiency  parameter 

P  =  Scattering  phase  function  normalized  to  unity  over  Att  steradians 
Tp  =  Temperature  of  the  plate 

fi  =  Angular  coordinate 

Because  z  <  <  L,  the  flat  plate  occupies  one  hemisphere,  and  the  angular  integral  in  Eq.(13) 
is  therefore  also  performed  only  over  a  hemisphere.  By  symmetry  and  the  normalization 
condition  on  P,  the  integration  of  P  over  a  hemisphere  equals  one-half.  The  radiance  scattered 
into  the  line  of  sight  is  therefore 

B(Tp)  (14) 
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where  Tp  is  the  temperature  of  the  flat  plate.  It  is  now  instructive  to  substitute  values 
appropriate  to  rocket  plumes.  We  use 


Qa  =  0.02, 

Q,  =  1.98, 

Td  =  2,300  K.  and 

Tp  =  1,300  K 


The  latter  is  a  representative  adiabatic  wall  temperature  at  the  nozzle  lip.  At  a  wavelength 
of  l.S  /un,  the  scattered  searchlight  emission  is  twice  the  directly  emitted  radiance. 

As  the  droplet  is  lifted  higher  above  the  plate,  the  searchlight  emission  decreases  for  two 
reasons.  First,  the  plate  no  longer  appears  infinite.  The  scattering  integral,  Eq.  (13),  is  limited 
to  a  smaller  solid  angle  which  approaches  the  asymptotic  limit  (L/z,  Ref.  2)  (i.e.,  inverse 
square  law)  for  large  z.  This  is  simply  the  radiative  exchange  view  factor  between  a  flat  plate 
and  a  distant  point.  Second,  the  scattering  angles  all  approach  a  limiting  value  of  90  deg. 
The  asymptotic  limit  of  the  searchlight  term  is  therefore 

11X2  Q,  P(90<’)  B(Tp)  (15) 

2.3  PLUME  RADIATION  MEASUREMENTS 

Selected  radiometer  and  calorimeter  data  from  two  motors  are  reviewed  to  illustrate  some 
of  the  rectirrent  features  of  exhaust  plumes  observed  at  the  AEDC  altitude  simulation  test 
chambers.  These  data  were  obtained  on  targets  of  opportunity  on  a  noninterference  basis 
during  routine  motor  quallHcation  and  testing.  They  are  not  a  complete  set,  and  at  no  time 
was  it  possible  to  perform  a  dedicated,  comprehensive  measurements  set  on  any  single  motor. 
However,  the  available  data  exhibit  characteristic  spatial,  temporal  ,and  wavelength  trends 
which  indicate  strong  searchlight  emission  at  the  nozzle  exit  plane. 

Even  for  small  motors,  the  exit  plane  radiance  does  not  scale  linearly  with  thrust.  Instead, 
it  tends  to  increase  more  slowly  than  thrust  after  ignition  and  then  to  continue  climbing  during 
the  late  bum,  somrtimes  even  after  the  thmst  is  level  or  tailing  off.  Figure  12  shows  the  radiance 
2  in.  downstream  of  the  nozzle  exit  plane  measured  by  a  narrow-ficld-of-view  calorimeter 
at  90-deg  aspect  angle  for  motor  A.  Although  the  thrust  peaks  at  35  sec,  the  radiation  does 
not  peak  until  65  sec.  In  addition,  the  initial  overshoot  in  thrust  near  10  sec  is  not  accompanied 
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by  a  corresponding  feature  in  the  radiance.  In  many  cases  there  are  no  further  data  from 
a  motor  firing,  and  the  cause  of  such  nonlinearity  remains  unknown.  For  this  particular  motor, 
however,  calorimeter  time  histories  were  also  measured  at  two  additional  downstream  locations 
(Fig.  13).  Although  all  three  calorimeters  were  aimed  at  the  centerline,  the  profiles  are  strikingly 
different.  In  particular,  as  one  progresses  downstream,  they  gradually  become  more  like  the 
thrust  curve.  This  is  consistent  with  a  localized  searchlight  enhancement  of  exit  plane  radiation, 
but  is  not  conclusive.  The  spikes  in  the  data  for  all  three  channels  are  attributed  to  the  repeated 
ejection  alumina  slag  from  the  motor  and  were  also  evident  in  the  video  as  bright  flashes. 

Spectrometer  data  provide  another  important  test  of  the  searchlight  mechanism.  In  the 
long  wavelength  infrared  (LWIR),  the  exhaust  particles  are  inefficient  scatterers  (Refs.  16 
and  17),  and  strong  searchlight  scattering  cannot  occur.  If  the  searchlight  hypothesis  is  correct, 
then  little  or  no  excess  radiation  should  appear  in  the  far  infrared.  Figure  14  confirms  this 
for  motor  B.  This  figure  was  generated  from  2-  to  14-/im  exit  plane  spectra  with  a  scan  time 
of  approximately  2  sec  per  spectrum.  It  compares  the  time  histories  of  the  short  wavelength 
infrared  (SWIR)  and  LWIR  emission  against  motor  thrust.  The  data  have  been  normalized 
at  the  second  time  point.  Whereas  the  LWIR  band  tracks  the  thrust  quite  closely,  the  SWIR 
band  deviates  from  it  as  it  did  for  motor  A.  This  result  indicates  that  the  anomalous  exit 
plane  radiance  profiles  are  associated  with  scattering  from  exhaust  AI2O3  in  the  plume.  In 
addition,  the  close  correspondence  of  the  LWIR  and  thrust  profiles  indicates  that  the  flow- 
field  and  particle  properties  remain  fmrly  steady  over  the  course  of  the  burn.  This  same 
measurement  was  unsuccessfully  attempted  on  a  larger  motor.  The  plume  was  opaque  and 
the  LWIR  radiance  profile  was  flat,  except  during  ignition  and  tailoff.  This  is  consistent 
with  the  high  LWIR  emissivity  of  AI2O3,  The  dual-wavelength  method  is,  therefore,  useful 
only  for  small  motors. 

Additional  information  was  sought  in  the  response  of  the  SWIR  radiance  to  fluctuations 
in  thrust.  The  fluctuations  under  consideration  are  fairly  slow  and  arise  from  normal  chamber 
pressure  excursions  as  the  area  of  the  burning  propellant  grain  changes.  The  bump  in  the 
thrust  profile  of  motor  A  near  10  sec  in  Fig.  12  is  a  good  example.  Since  the  plume  is  optically 
thin  with  respect  to  absorption,  the  direct  emission  from  the  plume  should  respond  almost 
instantly  (within  milliseconds)  and,  at  least  to  first  approximation,  linearly  to  thrust  changes. 
As  can  be  seen,  however,  this  is  not  the  case. 

To  better  quantify  this  observation,  the  thrust  and  radiance  data  of  Fig.  12  were  Fourier 
transformed  into  the  frequency  domain  (Fig.  15).  The  quantity  of  interest  is  not  the  particular 
shape  of  the  transforms,  but  rather  their  ratio,  as  this  gives  the  frequency  response  function 
of  the  SWIR  to  thrust.  If  all  of  the  radiance  were  due  to  direct  plume  emission,  the  response 
function  would  be  flat,  and  the  two  curves  would  overlay  each  other.  What  actually  occurs 
is  that  the  radiance  power  spectral  density  begins  to  fall  below  the  thrust  curve  above  roughly 
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0.2  Hz.  The  searchlight  hypothesis  affords  a  natural  explanation  for  this  in  terms  of  the  thennal 
lag  time  of  the  nozzle.  Specifically,  above  approximately  0.2  Hz  the  fluctuations  are  too  fast 
for  the  nozzle  temperature  to  track,  and  only  the  direct  plume  emission  contributes  to  the 
radiance  fluctuations.  Nozzle  thermal  lag  also  provides  a  reasonable  explanation  for  the  way 
in  which  the  SWIR  initially  lags  behind  the  thrust. 

In  the  case  of  motor  B,  a  calculated  nozzle  lip  radiance  time  history  B(t)  was  available 
(Fig.  IQ.  The  predicted  intensity  is  several  times  that  of  the  plume.  This  indicates  that  the 
efficiency  of  the  searchlight  scattering  mechanism  does  not  need  to  be  very  high  to  account 
for  the  data.  In  addition,  the  calculation  shows  that  the  thermal  lag  time  of  the  nozzle  lip 
is  consistent  with  the  observed  lag  times  of  exit  plane  radiance.  An  attempt  was  made  to 
see  how  well  the  observed  SWIR  time  history  could  be  explained  using  a  linear  sum  of  two 
functions,  one  proportional  to  the  thrust  (the  direct  plume  emission),  and  the  other  to  the 
calculated  nozzle  radiance  (the  searchlight  emission).  This  neglects  changes  in  the  searchlight 
scattering  properties  of  the  plume  with  thrust  and  so  is  only  approximate.  A  least-squares 
fit  to  the  two  components  gave  the  results  of  Fig.  17.  The  overall  profile  of  the  SWIR  data 
was  fairly  well  described  by  this  procedure.  However,  the  inferred  intensity  of  the  searchlight 
emission  was  more  than  twice  as  great  as  that  indicated  by  the  SWIR  versus  LWIR  comparison 
shown  previously  in  Fig.  14.  Since  the  nozzle  radiance  calculations  suffer  from  uncertainty 
in  the  high-temperature  properties  of  the  partially  pyrolyzed  carbon  phenolic  nozzle  material, 
the  SWIR  versus  LWIR  comparison  is  probably  more  reliable. 

3.4  MONTE  CARLO  ESTIMATi:  OF  SEARCHLIGHT  INTENSITY 

A  theoretical  basis  for  the  preceding  observations  can  be  established  using  the  Monte 
Carlo  method  for  radiative  transport.  This  method  provides  a  means  of  assessing  realistic 
plume  effects  such  as  multiple  scattering  and  inhomogeneity.  For  the  present  application, 
the  backward  Monte  Carlo  method,  sometimes  called  time-reversed  Monte  Carlo,  is 
computationally  much  more  efficient.  In  this  variant,  one  computes  absorption  rather  than 
emission.  By  Kirchoff  s  law,  the  two  are  equal  under  conditions  of  local  thermodynamic 
equilibrium.  Absorption  probabilities  are  estimated  by  ray  tracing  photons  with  small, 
randomized  step  size  through  the  plume.  The  probability  of  absorption  and  scattering  within 
each  step  is  determined  using  a  random  number  generator.  If  scattering  occurs,  a  new  random 
direction  is  chosen  consistent  with  scattering  phase  function.  The  ray  trace  is  continued  until 
the  photon  is  either  absorbed  or  escapes  from  the  plume  boundary.  Monte  Carlo  calculations 
for  the  emissivity  of  anisotropic  scattering  slabs  agreed  with  numerical  results  (Ref.  24)  to 
withm  the  expected  standard  deviation.  This  was  about  one  percent  for  10*  photons. 

The  Monte  Carlo  searchlight  calculation  requires  specification  of  the  spatial  distribution 
of  particles  vrithin  the  plume.  The  radial  distribution  of  particles  for  motor  A  was  based 
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upon  Abel  inversion  of  Nd:YAG  laser  transmission  measurements  performed  at  eight  radial 
locations  (Fig.  18).  The  quantity  determined  by  this  measurement  is  the  extinction  coefficient 
and  is  used  directly  by  the  Monte  Carlo  code.  Several  inversions  performed  approximately 
1  sec  apart  are  shown  for  two  time  segments,  one  near  the  beginning  and  one  near  the  end 
of  the  motor  burn.  Successive  inversions  are  shown  to  give  repeatable  results.  The  differences 
in  the  profiles  at  the  beginning  and  end  of  the  burn  are  attributed  to  differences  in  chamber 
pressure  and  the  internal  grain  configuration.  The  depletion  of  the  particle  extinction  at  the 
plume  centerline  is  not  due  to  spin,  as  the  motor  was  stationary.  The  highly  nonuniform 
nature  of  this  distribution  indicates  the  value  of  IR  camera  measurements.  The  particle 
distribution  at  other  axial  locations  was  scaled  from  th»e  results  assuming  a  conical  expansion. 

The  scattering  phase  function  was  approximated  by  the  Henyo^-Greenstein  function  (Refs. 
25  and  26).  This  is  an  analytical  approximation  to  the  Mie  scattering  phase  function.  It  lacks 
the  backscatter  lobe  near  180  deg,  but  it  is  otherwise  a  good  representation  of  the  full  Mie 
result,  especially  if  one  considers  a  distribution  of  droplet  sizes.  The  one  free  parameter  of 
this  function  specifies  the  degree  of  forward  scattering,  specifically  the  mean  cosine  of  the 
scattering  angle.  The  value  of  this  mean  angle  was  taken  to  be  0.5,  based  on  the  size  determined 
from  a  laser  transmissometer  and  Mie  calculations.  The  results  were  not  sensitive  to  the  precise 
value  of  this  parameter.  The  volume  scattering  albedo  and  scattering  phase  function  were 
assumed  constant  throughout  the  plume.  As  for  the  former,  the  actual  value  is  sufficiently 
close  to  unity  that  small  variations  do  not  affect  the  basic  scattering  behavior  of  the  plume. 
As  for  the  latter,  there  is  inadequate  information  to  make  a  more  detailed  choice. 

In  the  backward  Monte  Carlo  method,  the  searchlight  emission,  I,  is  given  by  the  product 
of  the  probability  of  absorption  on  the  nozzle  inner  wall,  dP/dz,  times  the  intensity  of  emission 
at  this  point,  B(z).  Because  the  nozzle  is  black  carbon  phenolic,  this  emission  was  taken  to 
be  the  Planck  function,  B,  at  the  wall  surface  temperature. 

I=jB(z)-^dz  (16) 

az 

The  quantity  dP/dz  can  also  be  regarded  as  the  differential  view  factor  between  the  radiometer 
and  the  inner  nozzle  wall.  The  region  of  integration  is  from  the  nozzle  exit  to  the  throat. 
If  the  nozzle  wall  temperature  is  assumed  constant  over  the  region  of  integration,  then  the 
Planck  function  can  be  taken  outside  the  integral.  The  searchlight  can  then  be  characterized 
by  a  view  factor,  F. 

dP 

F  =  J  dz  (17) 

dz 

Two  different  normalizations  for  view  factor  are  in  common  use.  The  one  used  here  is  EFi  = 

1. 
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The  probability  distributions,  dP/dz,  for  each  of  the  three  broadside  radiometers  used 
on  motor  A  were  determined  from  Monte  Carlo  histograms  of  100,000  photons  each  with 
a  bin  size  of  3  cm  in  the  axial  direction  (Fig.  19).  Each  distribution  contains  a  delta  function 
component,  due  to  directly  transmitted  photons,  and  a  continuous  bell-shaped  component 
due  to  scattered  photons.  The  directly  transmitted  component  increases  downstream  because 
plume  expansion  reduces  the  overall  scattering  losses.  Paths  terminating  to  the  left  of  the 
nozzle  exit  plane  (at  0  cm)  terminate  on  the  inner  nozzle  wall  and  therefore  contribute  to 
searchlight  emission.  The  integrals  of  the  curves  to  the  left  of  the  exit  plane  give  the  search' 
light  view  factors  [Eq.  (17)].  It  is  quite  high,  0.33,  at  the  2-in.  location.  This  is  2/3  of  the 
theoretical  maximum  value  of  O.S  at  the  exit  plane.  It  then  drops  rapidly  to  0.10  and  0.04 
at  the  17'  and  29'in.  stations,  respectively.  Comparison  with  the  radiometer  data  of  Fig.  13 
indicates  that  this  axial  failoff  is  consistent  with  the  data  trend.  As  shown,  most  of  the 
searchlight  emission  is  predicted  to  originate  primarily  from  the  lip  region  of  the  nozzle.  This 
is  consistent  with  the  fairly  slow  rise  time  of  the  excess  radiation  and  justifies  the  use  of  a 
net  view  factor  to  characterize  the  searchlight. 

The  implications  of  searchlight  emission  for  radiative  heating  were  briefly  alluded  to  in 
conjunction  with  Fig.  10,  which  indicated  that  the  forward  scattering  nature  of  the  plume 
particulates  would  tend  to  minimize  searchlight  emission  in  the  direction  of  the  vehicle.  This 
issue  is  pursued  further  in  Fig,  20,  which  shows  the  aspect  angle  dependence  of  the  search¬ 
light  view  factor  for  motor  A.  The  view  factors  are  computed  for  an  exit  plane  calorimeter 
aimed  at  the  plume  centerline.  The  method  and  flow-field  inputs  are  identical  to  those  used 
for  Fig.  19.  Between  90  deg,  where  data  were  taken,  and  30  deg,  which  is  more  relevant  to 
vehicle  heating,  the  searchlight  drops  by  a  factor  of  five.  Even  with  multiple  scattering,  the 
searchlight  does  not  contribute  appreciably  to  vehicle  heating.  Heating  rates  based  solely 
on  broadside  measurements  will  therefore  be  overly  conservative.  For  critical  applications, 
and  certainly  when  the  broadside  searchlight  is  strong,  a  more  realistic  heating  estimate  will 
result  if  the  searchlight  component  can  be  identified  and  removed  by  subtraction.  There  is 
a  second  and  perhaps  more  important  reason  for  subtracting  the  searchlight  term.  In  the 
current  industry  standard  method,  the  radiance  of  the  downstream  plume  is  extrapolated 
from  a  curve  flt  to  the  near-Held  (broadside)  data.  As  has  been  shown,  however,  the  spatial 
proHle  of  the  near-field  plume  radiance  can  be  confounded  by  searchlight  emission.  In 
particular,  if  the  searchlight  emission  is  strong,  then  a  curve  fit  to  the  near-field  radiance 
profile  may  have  only  limited  relevance  to  the  direct  plume  thermal  emission  which  one  actually 
needs  for  heat  transfer.  Since  the  searchlight  has  a  fast  axial  decay  rate,  moreover,  a  curve 
fit  made  without  allowances  for  searchlight  will  tend  to  underestimate  the  radiance  of  the 
far-field  plume. 
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3.5  SUMMARY 

The  prediction  of  in-fitght  radiative  heat  transfer  from  the  plumes  of  aluminized  solid- 
propellant  space  motors  is  difficult  because  typical  altitude  test  facilities  do  not  permit  direct 
measurements  of  the  heat  flux  at  the  desired  locations  and  viewing  angles.  Instead,  the  in¬ 
flight  heating  must  generally  be  extrapolated  from  a  limited  data  set  taken  near  the  nozzle 
exit  plane  and  basically  broadside  to  the  plume.  Measurements  of  exhaust  plume  radiation 
from  aluminized  solid-propellant  motors  at  the  AEDC  altitude  test  facilities  indicate  that 
both  the  intensity  and  spatial  distribution  of  this  radiation  are  confounded  by  searchlight 
emission.  The  broadside  exit  plane  radiance  does  not  scale  with  thrust,  even  for  small  motors. 
Instead,  it  lags  behind  the  thrust  curve  after  ignition  and  later  continues  to  rise  even  when 
the  thrust  is  level  or  tailing  off.  The  radiance  rise  time  is  several  seconds  and  is  comparable 
to  the  calculated  beating  rate  for  the  nozzle  lip.  At  downstream  locations,  on  the  other  hand, 
the  radiance  tends  to  track  the  thrust  curve  more  closely.  The  anomalous  exit  plane  radiance 
profile  occurs  only  in  the  SWIR  where  alumina  is  a  strong  scatterer  and  is  indiscernible  in 
the  LWIR  where  it  is  a  poor  scatterer.  Finally,  examination  of  the  radiance  and  thrust  power 
spectral  densities  in  the  frequency  domain  indicates  the  presence  of  two  distinct  components, 
one  which  responds  rapidly  to  thrust  fluctuations,  and  the  other  which  responds  with  a  time 
constant  of  several  seconds. 

These  observations  are  attributed  to  searchlight  emission.  Specifically,  they  are  attributed 
to  hot  nozzle  wall  radiation  which  is  scattered  into  the  broadside  direction  by  the  exhaust 
alumina  particulates.  The  strength  and  the  axial  falloff  of  the  searchlight  emission  are 
supported  by  Monte  Carlo  radiative  transfer  calculations  using  the  measured  plume  extinction 
profile.  It  is  concluded  that  searchlight  emission  is  much  stronger  than  has  been  previously 
supposed.  Although  it  contributes  strongly  to  the  broadside  exit  plane  radiance,  it  is  probably 
not  important  at  the  smaller  aspect  angl»  relevant  to  vehicle  heating.  A  more  realistic  estimate 
of  the  near-field  plume  radiance  will  result  if  the  searchlight  component  can  be  identified 
and  subtracted  from  the  broadside  radiance  data.  In  addition,  an  improved  extrapolation 
of  the  far -field  plume  radiance  will  also  result,  allowing  the  industry  standard  method  to 
produce  more  accurate  estimates  of  heat  flux.  If  these  corrections  can  be  performed,  the 
current  industry  standard  method  should  then  provide  accurate  estimates  of  the  in-flight 
radiative  heating  load.  To  this  end,  several  techniques  for  the  identification  and  quantification 
of  searchlight  emission  were  attempted.  Of  the  methods  tried  thus  far,  the  best  indicator 
is  probably  comparison  of  the  exit  plane  versus  downstream  radiance  time  histories. 
Comparison  of  the  SWIR  and  LWIR  time  histories  appears  to  be  a  useful  alternative  for 
small  motors,  for  which  the  LWIR  alumina  absorption  is  optically  thin.  It  has  the  additional 
advantage  of  being  applicable  even  when  only  a  small  length  of  plume  is  visible.  The  increased 
opacity  of  large  motor  plumes  prevents  its  application  for  these  motors. 
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4.0  SUMMARY 

The  analyses  of  the  AlCl  absorption  feature  and  the  searchlight  effect  observed  in  solid- 
propellant  rocket  plumes  during  chamber  tests  at  simulated  altitudes  show  that  these 
phenomena  are  of  particular  interest  in  the  characterization  of  plumes.  The  AlCl  absorption 
analysis  indicates  that  this  feature  may  be  of  use  for  plume  diagnostics,  although  additional 
data  are  required  before  a  better  fit  between  calculated  and  measured  spectra  can  be  obtained 
(Ref.  27).  The  searchlight  analysis  indicates  that  searchlight  emission  is  a  strong  component 
in  the  near-fleld  plume  and  needs  to  be  taken  into  account  when  attempting  to  extrapolate 
far-fleld  plume  radiance  from  near-field  measurements  (Ref.  28). 
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Figure  1.  Schematic  of  soUd-propellant  rocket  motor  test  configuration  used 
for  UV  spectral  measurements. 


27 


AEDC-TR-92-1 2 


Figiire  2.  Modaste-resolution  mld-TJV  qiectra  firom  Mdtd-propdlant 
rocket  motor  plume. 


Figure  3.  High-resolution  UV  spectra  from  solid-propellant  rocket 
motor  plume. 
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b.  Pressure  profile 

Figure  4.  Predicted  radial  profiles  for  solid-propellant  rocket  motor 
at  axial  location  where  UV  spectra  were  measured  (stan¬ 
dard  plume  flow-field  code  using  three  particle  groups). 


29 


AEDC-TR-92-12 


PLUME  CROSS  SEQION 


Figure  5.  Physical  model  of  solid-propellant  rocket  motor  plume 
UV  spectral  measurement  experiment. 


Figure  6.  Comparison  of  computed  with  measured  spectra  of  SRM 
plume  radiation  at  predicted  plume  conditions. 
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Figure  7.  Comparison  of  compated  with  measnred  spectra  of  SRM 
plume  radiation  at  “best  fit”  conditions. 


Q  RADIOMETER 


Figure  8.  Sketch  of  searchlight  emission  where  small  angle 
scattering  hy  plume  AI203  droplets  and  particles 
redirects  nozzle  emission  toward  broadside 
radiometers. 
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SCAHERING  ANGLE,  DEG 


Figure  10.  Scattering  phase  function  of  AI2O3  droplet  at  three 
wavelengths. 
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Figure  11.  Geometry  of  gedanken  experiment  showing  radiance  of 
a  reflecting  and  emitting  sphere  suspended  above  a  hot 
plate. 


Figure  12.  Motor  A  thrust  and  SWIR  radiance  time  histories. 
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Figure  IS.  Power  spectral  density  comparison  of  thrust  versus 
short  wavelength  infrared  (SWIR)  radiance. 


Figure  16.  Calculated  radiance  of  motor  B  nozzie  iip  at  two 
wavelengths. 
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Figure  17.  Decomposition  of  motor  B  radiance  into  searchlight  and 
plume  emission  based  upon  calculated  nozzle  lip 
temperature  time  history. 


Figure  18.  Extinction  coeffident  for  motor  A  where  inversions  at 
approximately  1  sec  intervals  are  shown  near  the  begin> 
ning  and  the  end  of  the  bum. 
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-100  -60  -20  20  60  100  140 


AXIAL  POSITION,  cm 

Figure  19.  Monte  Carlo  calculation  of  searchlight  for  motor  A  for 
radiometers  located  at  2-,  17-,  and  29-in.  axial  position 
where  the  nozzle  exit  plane  is  at  0  cm  and  the  table  gives 
searchlight  view  factor  for  each  radiometer. 


ASPEQ  ANGLE 

Figure  20.  Monte  Carlo  calculation  of  searchlight  view  factor 
versns  viewing  aspect  angle  where  radiometer  loca¬ 
tion  is  2  in.  downstream  of  nozzle  exit  plane. 
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Table  1.  Bandheads  of  the  AlCl  A^Il  -  System 


X 

1 

V,V 

1 

V,V 

Z708.9  R 

2 

R  7,9 

2625.0  R 

5 

Q3,3 

2702.3  R 

3 

R  6,8 

2622.4  R 

4 

R3,3 

2696.4  R 

3 

R  5,7 

2620.0  R 

4 

Q  2,2 

2692.8  R 

5 

Q4,6 

2618.2  R 

3 

R  2,2 

2685.7  V 

6 

Q  2,4 

2617.0  R 

4 

Ql,l 

2683.1  V 

6 

Q  1,3 

2614.4  R 

8 

Q0,0  R1,1 

2681. 1  V 

4 

Q  0,2 

2610.2  R 

6 

R  0,0 

2649.7  V 

4 

Q  1,2 

2606.7  R 

2 

R6,5 

2647.5  V 

6 

Q  0,1 

2600.7  R 

3 

R5,4 

2644.9  R 

2 

R  7,7 

2695.4  R 

2 

R  4,3 

2638.1  X 

3 

R  6,6 

2590.8  R 

2 

R  3,2 

2632.8  R 

3 

0  5,5 

2586.7  R 

2 

R2,1 

2632.2  R 

3 

R5,5 

2564.3  R 

1 

R4,2 

2627.8  R 

4 

0  4,4 

2559.6  R 

1 

R  3,1 

2627.0  R 

3 

R4,4 

2555.5  R 

1 

R  2,0 

NOTE:  THE  LETTERS  R  AND  V  AFTER  THE  WAVELENGTH  INDICATE 
THE  DIREQION  OF  DEGRADATION  OF  THE  BAND,  WHEREAS 
THE  NATURE  OF  THE  HEAD  (R  AND  Q)  IS  INDICATED  BEFORE 
THE  VIBRATIONAL  QUANTUM  NUMBERS. 

REF.:  THE  IDENTIFIUTION  OF  MOLECUIAH  SFEaU,  R.  W.  B. 
PEARSE  AND  A.  G.  GAYDEN,  CHAPMAN  AND  HAU  LTD.,  1950, 
p.  47. 
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Table  2.  Molecular  Constauts  for  the  Ground  and  Resonant  States  of 
AI3SC1  (cm-i) 


x’i:+ 

A'n 

Te 

0 

38254.0 

481.30 

449.96 

1.95 

4.37 

— 

-0.216 

Be 

0.24393 

0.259 

“e 

O.ODIAlIt 

0.006 

Da 

2.503  X  10-^ 

— 

-0.0053  X  10-7 

... 

BASED  UPON  Q-BRANCH  BANDHEADS 

T  +4.697  X  10-‘(v  +  '/j)2  -  4.7  x  lO"’  (v  + 


Table  3.  Einstdn  Coefficients  for  the  Ain  Transition  in  AlCi 

(Ref.  7). 


.\y" 

0 

1 

2 

3 

4 

5 

6 

7 

8 

0 

1933 

— 

—  - 

... 

... 

... 

... 

... 

... 

1 

... 

1890 

— 

— 

... 

... 

... 

... 

2 

— 

7.92 

18.32 

... 

... 

... 

... 

... 

3 

— 

— 

3.3 

1745 

12 

16.6 

— 

... 

... 

4 

— 

— 

— 

92 

1578 

37 

31 

... 

... 

5 

... 

... 

... 

... 

200 

1344 

77 

54 

... 

6 

... 

... 

— 

— 

— 

366 

1028 

121 

87 

7 

... 

— 

... 

... 

— 

— 

574 

658 

149 

8 

... 

... 

... 

— 

— 

— 

18 

766 

301 

ALL  VALUES  SHOULD  BE  MULTIPLIED  BY  lO^  set-l  (REF.  6} 

VALUES  NOT  GIVEN  ARE  LESS  THAN  0.5  PERCENT  OF  VALUE  AT  (0,0) 
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Table  4.  Constants  for  Dunham  Level  Expression  of  Energy  Level  of 
the  X'L  State  of  AlCl 


Moletular  Constants  for  AICI(Xli;]  used  In  Term  Energy  Expression: 

T,J  “  plmlv  +  Wj'Jmy  +  1)" 

l,in 

^Im 

1 

m 

AI35CI 

Ai37a 

O 

D 

O.OOOOODOOOOOOE+00 

O.OOOOOOOOOOOOE+00 

n 

n 

2.439300631 168E~0I 

2.381908103047E-01 

0 

2 

-2.5016234841311-07 

-2.3846883S5076E-07 

0 

3 

-4.2478467I78I6E-14 

-5.6318344541 37E-14 

0 

4 

9.D82196640775E-19 

8.257101917000E-19 

a 

D 

4.817765661053E+02 

4.760760032685E+02 

D 

n 

-1.611082122620E-03 

~1. 554551 027648E-03 

1 

2 

4.398230380075E-10 

4.161398270080E-10 

1 

3 

O.OOOOOOOOOODOE+00 

O.OOOOOOOODOOOE+00 

2 

O 

-2.10181 7969356E+00 

-2.052568958832E+00 

2 

1 

4.6920560 13898E-06 

4.47253471 8003E-06 

2 

2 

7.160332258952E-12 

8.149446348886E-12 

3 

D 

6.6381 64756386E-03 

6.42086760831 9E-03 

3. 

-5.291 195847775E-05 

-4.98542057BO00E-O9 

D 

D 

-2.023B33735279E-05 

-1.929715135000E-fl5 

VALUES  ARE  IN  UNITS  OF  CM-1 
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